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STRAIN-HARDENING  AND  INTERACTION  EFFECTS  ON 
THE  GROWTH  OF  VOIDS  IN  DUCTILE  FRACTURE 

ABSTRACT 

The  growth  and  coalescence  of  voids  is  a  common  mechanism  of  fracture  in 
ductile  materials.  Analytical  work  on  the  problem  to  date  has  dealt  mainly  with 
isolated  voids  in  perfectly  plastic  materials,  so  that  strain  hardening  and 
interactions  between  neighboring  voids  have  been  neglected.  These  features  of 
void  growth  are  examined  here,  but  only  for  a  simple  geometrical  configuration. 

In  particular,  the  growth  of  single  infinitely  long  cylindrical  voids  in  bodies 
of  rigid-plastic,  strain-hardening  material  is  considered.  Bodies  with  both 
finite  and  infinite  dimensions  normal  to  the  void  surface  are  included  in  the 
analysis.  The  exact  relation  among  the  pertinent  variables:  transverse  stress, 
axial  strain,  hardening  coefficient,  void  strain  and  void  growth  rate  is  pre¬ 
sented.  Solution  via  a  bounding  technique  is  given  for  two  general  cases.  The 
first  case  is  that  of  an  imposed  constant  void  growth  rate  and  the  second  case  is 
an  imposed  constant  transverse  stress.  The  results  show  a  decelerating  effect 
of  hardening  on  void  growth.  Application  to  the  ductile  fracture  problem  of 
void  growth  in  the  neck  of  a  tensile  specimen  demonstrates  the  accelerating 
effect  of  void  growth  due  to  interactions  between  voids. 


Introduction 


Ductile  fracture  as  a  consequence  of  the  growth  and  coalescence  of  voids 
has  been  observed  by  many  experimenters  including  Rogers  [1]  ,  Gurland  and 
Plateau  [2],  and  Bluhm  and  Morrissey  [3].  Voids  nucleate  by  the  cracking  of 
inclusions,  the  decohesion  of  inclusions  from  the  matrix,  and  reportedly  also 
in  the  matrix  metal  near  the  inclusions  [2].  McClintock  and  O' Day  [4]  discount 
the  importance  of  the  last  nucleation  site  in  their  finding  that  dislocations 
have  a  small  effect  on  fracture  strain  of  the  bubble-raft  analogue  of  a  crystal. 
Once  the  voids  appear  they  grow  and  coalesce  by  the  plastic  deformation  of  the 
matrix  which  in  this  paper  is  considered  a  problem  in  continuum  plasticity. 
McClintock  [5]  justifies  this  approach  as  reasonable  when  the  inclusions  in  the 
material  are  large  enough  to  be  seen  under  an  optical  microscope. 

The  fracture  of  a  simple  tension  specimen  commences  with  the  nucleation 
of  fine  cracks  and  voids  in  the  necked  region  along  the  axis.  Rogers  [1]  states 
that  the  flat  portion  of  the  cup  and  cone  fracture  is  a  direct  consequence  of  the 
coalescence  of  voids  which  have  growth  under  the  influence  of  the  prevailing  tri- 
axial  stress  system.  The  voids  are  most  dense  at  the  axis  of  the  neck  where 
the  triaxiality  of  stress  is  greatest  indicating  a  direct  relation  between  void 
nucleation  and  growth  and  triaxiality. 

The  analytical  approach  to  the  problem  is  not  obvious  because  the  mechanism 
of  coalescence  is  not  understood.  A  first  attempt  would  ignore  the  interaction 
between  voids  and  consider  each  void  as  growing  in  an  infinite  domain  until  the 
radius  has  growth  to  one-half  the  spacing  between  voids.  McClintock  [6j  used 


The  numbers  in  square  brackets  refer  to  the  list  of  references  at  the  end 
of  this  report. 


this  approach  in  his  void  growth  work.  The  approach  in  this  paper  again 
ignores  interaction  as  such  but  considers  each  void  as  in  its  own  secluded 
domain  of  transverse  dimension  equal  to  the  average  spacing  between  voids. 

At  some  small  void  density  it  is  valid  to  consider  the  void  domains  as  in¬ 
finite  ;  as  growth  proceeds  the  effects  of  the  void  growth  will  be  "felt"  in 
the  neighboring  voids'  domains  so  that  finite  body  analysis  begins  as  an 
approximate  method  of  describing  the  acceleration  in  growth  due  to  the  interaction. 

There  has  been  no  published  work  known  to  the  writer  which  solves  exactly 
the  problem  of  interaction  of  voids  (however,  McClintock  [5]  presented  an 
approximate  analysis  of  coalescence  by  shear  bands  between  voids);  in  fact, 
there  has  not  been  any  work  reported  on  the  exact  solution  of  non-interacting 
voids  in  strain-hardening  materials.  McClintock  [6]  presented  a  criterion  of 
ductile  fracture  by  the  growth  of  non-interacting  voids  in  infinite  domains  of 
rigid  strain-hardening  material.  The  void  growth  behavior  for  the  strain- 
hardening  material  is  not  exact  however;  it  is  obtained  from  an  extrapolation 
between  the  behavior  in  perfect  plastic  and  viscous  materials .  Rice  and 
Tracey  [7]  presented  a  variational  treatment  of  void  growth  in  infinite  domains 
of  any  rigid  plastic  material  but  only  perfect  plastic  applications  were  pre¬ 
sented.  Hence  an  exact  strain-hardening  treatment  is  desirable  and  here 
presented. 

Formulation  of  Problem 

In  a  previous  paper  [7]  the  growth  of  isolated  cylindrical  and  spherical 
voids  in  infinite  bodies  of  perfectly  plastic  material  was  treated.  Here  the 
growth  of  a  cylindrical  void  in  "finite"  (infinite  length  and  finite  transverse 
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dimensions)  and  infinite  bodies  of  rigid-plastic,  strain-hardening  material  is 
considered.  Fig.  1.  The  bodies  are  cylindrical  in  shape  with  an  axis  of  in¬ 
finite  length  and  a  radius  b,  b  having  any  value,  finite  or  infinite.  A 
cylindrical  void  of  radius  rQ  is  oriented  along  the  z  direction  of  the 
(r,6,z)  polar  coordinate  system  defined  in  Fig.  1.  A  strain  is  imposed  on  the 
body  in  the  axial  direction  at  a  rate  e  while  the  radial  stress  (o^)^  is 
applied  at  the  external  radius  b  .  The  axial  direction  is  the  direction  of 
maximum  principal  stress  and  strain-increment  in  all  cases  considered.  The 
problem  lies  in  finding  the  expansion  of  the  void  as  a  function  of  the  imposed 
stress  and  strain. 


The  incompressibility  condition  gives  the  velocity  field  [7] 
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The  void  growth  rate  r0^r0  must  he  known  to  completely  specify  the  velocity 
field.  Note  that  at  large  values  of  r  the  velocity  at  r  is  not  influenced 
by  the  value  of  r0/r0  •  If  the  body  has  a  large  enough  radial  dimension  so 
that  the  void  growth  is  not  "felt"  at  some  large  value  of  r  within  the  body, 
then  the  body  can  be  considered  infinite. 


The  relation  among  void  growth  rate,  imposed  stress  and  strain,  and  strain 
history  is  formulated  through  the  equation  of  equilibrium 
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and  the  stress-incremental  strain  relation  which  for  an  isotropic  Mises 
material  is 
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where  s..  is  the  deviatoric  stress  tensor,  e..  is  the  strain-rate  tensor 
13  13 

and  t  is  the  hardening  parameter  which  is  a  function  of  the  deformation  his¬ 
tory  and  equal  to  the  equivalent  shear  flow  stress.  Obtaining  e...  from 
eq.  (1),  substituting  into  eq.  (2),  dividing  both  sides  by  (the  equivalent 

shear  flow  stress  at  the  boundary),  and  integrating  over  the  body  gives 
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A  power  law  relation  between  t  and  the  generalized  strain  y  is 
assumed.  In  Lagrangian  coordinates  (radial  coordinate  denoted  by  R)  y  can 
be  written  as 


y(e,R) 


e 

y(n ,R)dn 
o 


(5) 


where 


The  strain-rate  field  e ^  is  found  from  the  velocity  field  eq.  (1)  in 
Eulerian  coordinates.  The  relation  between  Eulerian  (r)  and  Lagrangian  (R) 
coordinates  is  found  through  incompressibility: 

(r/r  )2  =  1+e'j  (1+2  ~)dn  [(R/R  )2  -ll  .  (6) 

o  ;o  o  L  °  J 

Using  eq.  (6)  to  obtain  y(R)  and  subsequently  transforming  back  to  Eulerian 
coordinates  to  obtain  y(r)  and  substituting  into  eq.  (4)  with  the  notation 
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gives  the  final  form  of  the  integral  equation  in  e,  rQ/b,  £ 
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Bounding  the  Flow  Stress 


The  flow  stress  t  appearing  in  eq .  (7)  is  a  cumbersome  quantity.  It  is 
a  function  of  position,  void  strain,  imposed  strain,  and  the  complete  strain 
history : 


In  an  effort  to  simplify  the  problem,  bounds  to  x  were  derived  by  obtaining 
bounds  to  y  which  hereafter  will  be  called  ya  .  y3  can  be  written  in 
Eulerian  coordinates  as 
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It  can  be  thought  of  as  the  arclength  of  the  actual  strain  path  in  e ,  H 
space  just  as 
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is  the  arclength  of  a  path  C  in  X^?  space, 
which  have  a  slope  of  constant  sign. 
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For  monotonic  paths ,  paths 
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This  point  is  illustrated  in  Fig.  2.  The  upper  bound  is  a  square  path  while 
the  lower  bound  is  a  straight  line  path.  The  same  bounds  can  be  applied  to  the 
arclength  of  the  actual  straining  path  ya  as  long  as  the  straining  process 
gives  a  monotonic  path  (of  course  the  path  can  always  be  divided  into  piecewise 
monotonic  portions ) : 


YU  =  |/3  e|  +  |*n((e"H(e)  -l)x2+l)|  }  ya  *  /3e2+(*n{(e‘H(E^  -l)x2+l})2  =  yl 

(12) 

The  bounds  are  path  independent;  this  is  a  very  significant  simplification 
from  eq.  (8).  They  were  employed  in  the  following  with  a  large  saving  in  com¬ 
putation  time. 


Solution  of  Integral  Equation  for  an  Infinite  Body  with  a  Constant  Void 
Growth  Rate 

When  the  void  growth  rate  is  constant,  H  =  constant,  the  integral 
eq.  (7)  has  its  simplest  form.  Although  the  physical  situations  in  which  a 
constant  growth  rate  is  encountered  are  not  numerous ,  solution  of  the  equation 
for  this  case  gives  a  feeling  for  the  void  growth  behavior  with  the  minimum  of 
computational  effort.  For  convenience,  the  body  containing  the  void  will  be 
chosen  to  have  an  infinite  radius.  The  resulting  integral  equation  is 
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FIGURE  2  BOUNDS  TO  MONOTONIC  STRAINING  PATHS 
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For  a  given  constant  value  of  H  eq .  (13)  is  solved  for  (o^/x^e)  ,  the 
stress  state  at  infinity  necessary  to  maintain  H  .  For  this  problem  the 

cl  £ 

bounds  to  y  can  be  applied  to  obtain  bounds  to  (o^/tJ^e)  .  By  using  y 
a  lower  bound  to  (o^/x)^  is  obtained;  by  using  yU  an  upper  bound  to 
(o  / t)  is  obtained: 
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The  quantity  H(e)  is  known  in  this  problem  to  equal  He  .  The  cumbersome 
double  integral  (13)  has  been  replaced  by  two  single  integrals.  The  computation 
time  of  obtaining  the  bounds  to  (o  /t)  (e)  is  many  times  smaller  than  the  time 

^  CO  v 

necessary  to  obtain  the  actual  (o^/tMe)  .  Fig.  3  shows  how  the  actual 
curve  (found  from  solution  of  eq .  (13)  falls  between  the  bounds  for  the  constant 
growth  rates  of  fQ/ro  =  5  and  lOe  . 
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Plots  of  (0^/t)^’  are  given  in  Fig.  4  for  growth  rates  of  5,  10,  20 
and  30e  for  N  =  .1  .  Fig.  3  is  the  same  for  N  =  .3  .  The  curves  were 
calculated  to  values  of  e  which  would  cause  a  geometric  strain  of  the  void 
equal  to  fifty. 

It  is  interesting  that  the  stress  state  required  to  maintain  the  imposed 
growth  rate  asymptotically  decreases  to  the  value  required  if  the  material  were 
perfectly  plastic  (N  =  0).  The  deviation  from  the  perfectly  plastic  stress 
state  at  a  value  of  e  decreases  with  value  of  imposed  growth  rate  and  hardening 
coefficient.  This  suggests  that  for  situations  which  induce  small  growth  rates, 
in  the  order  of  e  ,  perfect  plasticity  would  be  an  acceptable  idealization. 

This  point  is  demonstrated  in  the  last  section.  Another  implication  of  the 
curves  of  Figs.  3-4  is  that  if  the  stress  state  were  imposed  constant,  the 
void  growth  rate  would  continually  increase  to  the  perfectly  plastic  growth  rate. 
The  next  section  investigates  this  possibility. 


Solution  of  Integral  Equation  for  an  Infinite  Body  with  a  Constant  Imposed 
Stress  State 


When  the  imposed  stress  state  is  constant  the  integral  eq.  (7)  has 
the  form 
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where  k  is  the  imposed  stress  (c^/t^  •  The  yield  stress  r  depends  upon 

the  imposed  strain  e  ,  position  x  ,  and  the  entire  history  of  the  void  strain 
which  is  represented  by  H(n)  .  The  problem  is  to  obtain  the  function  H(n)  or 
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FIG. 4  BOUNDS  TO  (G> /T)®  [€]  FOR  IMPOSED  GROWTH  RATES 
OF  5,  10,20,30.  N  =  .1 


13 


just  as  well  H(n)  which  gives  H(n)  upon  integration.  Again  bounds  to  i 
will  be  employed. 

The  proportional  strain  and  "square"  strain  bounds  to  the  actual  strain 
path  will  be  used  to  release  t  of  its  path  dependence: 

T*(e,  H(e),  x)  *  x(e,  H(n),  x)  $  xu(e,  H(e),  x)  (16) 


£  u  jj, 

where  t  and  x  are  from  eq.  (12).  The  use  of  t  would  render  an  upper 

bound  to  H(e)  while  tu  would  render  a  lower  bound  to  H(e)  if  x11’8,  could 

be  evaluated.  They  cannot  be  evaluated  since  they  depend  upon  H(e),  the  intrinsic 

.  .  •  £,  u 

unknown  of  the  problem.  Since  the  derivatives  of  x  and  t  with  respect  to 
H(e)  are  positive,  bounds  to  H(e)  can  still  be  obtained  by  replacing  H(e) 

l  £  u 

in  t  with  some  lower  bound,  H  (e),  and  replacing  H(e)  in  t  with  some 
upper  bound,  Hu(e)  .  The  resulting  bounds  on  H(e)  ,  Hu  and  H2,  ,  will  bound 

U  £ 

the  bounds  that  would  result  from  the  use  of  H(e)  in  x  ,  t  .  Obvious  can¬ 
didates  for  HU(e)  and  H2'(e)  are 
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It  should  be  noted  that  any  bounds  to  H(e)  can  be  employed. 


Substituting  into  the  integral  equation  (13)  yields  two  single  integral 
equations : 
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Equations  (18)  are  two  equations  in  the  unknowns  H(e)  *  .  They  are  inter- 

related.  The  first  yields  H  (e)  but  only  after  knowing  H  (e)  which  is  found 

*  Z 

from  the  second  equation  by  integrating  H  (n)  .  Similarly  the  second  equation 
requires  H  (n)  before  yielding  H  (e)  . 

At  e  =  0  the  exact  equation,  eq.  (15),  reduces  to 


*  ck  JL  \1 

Solving  for  H  (0)  and  using  it  to  evaluate  t  (0)  and  t  (0)  allows  solution 

*  Vi  •  £ 

of  eq.  (18)  for  H(0)  and  H(0)  .  Thereupon  the  bounds  are  iteratively 

evaluated  at  small  increments  of  e  .  In  the  plots  presented  the  increment  in 

e  was  adjusted  with  k  so  as  to  instigate  a  void  true  strain  increment  of  .05. 

Plots  of  the  bounds  to  r  /er  were  obtained  for  (o  /t)  values  of  3,  3.5,  4, 

o  o  r  co 

and  5.5  for  hardening  coefficients  of  .1  and  .3  ,  Figs.  5-7.  Actual 
solution  of  eq.  (15)  was  obtained  for  imposed  stress  states  of  3  and  5.5  , 
Figs.  6,7. 

The  spread  between  bounds  is  greatest  for  high  stress  states.  For  small 
imposed  stresses  like  those  encountered  in  the  neck  of  a  tensile  specimen  the 
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bounds  are  very  close.  In  all  cases  the  void  growth  rate  increases  asympto¬ 
tically  to  the  perfect  plastic  growth  rate  as  the  longitudinal  strain  z 
increases.  The  degree  to  which  strain-hardening  decelerates  void  growth  is 
directly  related  to  the  degree  of  triaxiality  and  the  value  of  the  hardening 
coefficient . 

The  growth  rate  predicted  from  McClintock’s  approximate  work  [6]  is  given 
in  a  few  cases  on  Figs.  5,7.  His  equation  (27)  for  this  case  of  equal  trans¬ 
verse  stress,  put  in  present  notation,  relates  void  growth  rate  and  stress  as 

r  /zv  =  {(l-n)(o  /t)  >  -l] 

o  o  2  1-n  r  ® 


The  transient  effect  of  void  growth  rate  with  strain  z  is  not  included 
in  this  approximation,  a  very  important  feature  as  seen  in  the  Figs.  5-7. 


When  the  strain-rate  along  the  axis  of  the  body  of  Fig.  1  is  reversed 
in  direction  so  that  the  z-direction  becomes  the  direction  of  minimum  prin¬ 
cipal  stress  and  strain-rate,  the  only  difference  in  result  from  the  present 
case  is  that  the  growth  ratio  is  faster  by  the  amount  z  .  The  integral 
equation  is  the  same  for  both  cases.  The  only  difference  is  in  the  relation 

between  H  and  r  /r  : 

o  o 


r  /r 
o 


o 


H  -t  e 
2 


The  minus  sign  holds  for  the  present  case,  the  positive  sign  when 


z 


is  reversed. 
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Void  Growth  in  the  Neck  of  a  Tensile  Specimen 


The  necking  of  a  tensile  specimen  introduces  triaxiality  of  stress  in  the 
necked  region.  Bridgman  [8]  found  that  the  strain  was  very  nearly  uniform  at 
the  necked  section,  so  that  the  state  of  stress  in  the  neck  consists  of  a  uni¬ 
axial  stress  equal  to  the  flow  stress,  vZ  t  ,  plus  a  hydrostatic  tension, 
a  .  The  latter  is  given  approximately  in  terms  of  the  neck  radius,  a  ,  radius 
of  curvature  of  the  neck  surface ,  R  ,  and  distance  from  the  axis ,  r  ,  as 

a  =  (/3  t)  In  (1  +  a/2R  -  r2/2aR)  .  (20) 

In  an  attempt  to  predict  the  ductility  of  tensile  specimens  for  various  ini¬ 
tial  void  densities  the  interior  of  the  necked  region  was  assumed  to  contain  a 
number  of  cylindrical  voids  of  radius  r^  and  spacing  2b  ,  with  axes  aligned 
with  the  tension  direction.  The  initial  void  diameter  to  spacing  ratio  was  taken 
as  the  cube  root  of  the  volume  fraction  of  voids  present  in  the  material  under  con¬ 
sideration.  Each  void  is  assumed  to  grow  independently  in  its  own  finite  domain 
of  radius  b  .  When  the  void  radius  grows  to  a  value  equal  to  b  there  is  coal¬ 
escence  and  internal  fracture  is  complete;  incipient  shear,  in  the  scheme  of  Bluhm 
and  Morrissey  [3],  follows  to  finalize  the  fracture  with  a  comparably  small 
additional  strain. 

Since  voids  are  seen  to  be  most  dense  at  the  axis  of  the  necked  region  [3], 
the  stress  state  there  will  be  used  as  that  stress  state  imposed  on  the  cell 
walls.  To  obtain  a  relation  between  a/R  and  the  longitudinal  true  strain  e 
the  a/R  vs  £  plot  of  Bridgman  [8]  was  utilized.  The  data  was  collected  using 
steel,  bronze,  and  copper  tensile  specimens .  The  spread  in  data  is  considerable 
(as  would  be  expected  since  the  necking  process  varies  with  material)  but  for  a 
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first  attempt  the  relation  obtained  from  fitting  a  curve  to  the  data  was  used: 


a/R  =  2.93  [1-  exp  {-(e  -  .l)/3}] 

for  e  >  .  1 


(21) 


With  this  the  imposed  radial  stress  can  be  written  as  a  function  of  e  : 


°r(e)  =  (^  t)  in  [1  +  —  ?3  (1-  exp  [-  (e  -  ,l)/3]}] 

for  e  >  .1 


(22) 


The  quantity  ( /3  t)  is  the  uniaxial  flow  stress  present  in  the  neck  along  with 
the  hydrostatic  stress.  Since  t  is  the  flow  stress  in  shear,  which  is  constant 
across  the  neck  (8),  it  was  evaluated  via  the  power-law  relationship  from  the 
strain  state  far  from  the  void  conglomeration  at  the  axis: 


t  =  t  «  (/3  e)' 

05 


(23) 


The  integral  equation  (7)  takes  the  form 


/3  in 


r  2.93  ,,  r  , 

1  +  — (1-  exp  [-  (e 


•  l)/3]} 


2 

rl 

Y 

r 

/3  e 

H(e)  xdx 
/3e2+x4H2(e) 


(29) 
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The  void  spacing  is  changed  only  by  the  nominal  transverse  strain 


b(e)  =  b(0 )  exp  {-c/2} 


so  that  the  radius  ratio  is  given  by 


exp  {  ~ 


H(n)dn) 


(25) 


The  hardening  coefficient  of  a  material  is  influenced  by  the  density  of  inclu¬ 
sions,  particles  and  voids.  The  data  of  Edelson  and  Baldwin  [91  relating  N  and 
the  volume  portion  of  second  phase  particles  (including  voids)  of  two-phase 
copper-base  alloys  was  used  in  this  analysis.  The  bounds  to  y  ,  eq.  (12),  can  be 
utilized  with  the  modification  of  using  bounds  to  H(e)  as  in  the  previous  sec¬ 
tion.  The  equations  were  solved  for  initial  void  volume  fractions  of  0.02  to 
0.2  .  The  fracture  strain,  ,  at  which  the  voids  coalesce  is  that  strain  when 
rQ/b  =  1  .  In  practice  at  some  radius  ratio  less  than  1  the  growth  rate  takes  a 
sudden  surge  which  implies  coalescence  at  the  next  increment  of  strain.  The  radius 
ratio  just  before  coalescence  varied  between  2/3  and  3/4  .  The  resulting 
ductility-volume  fraction  curve  is  plotted  in  Fig.  8  along  with  the  experimental 
curve  of  Edelson  and  Baldwin  [9],  the  curve  offered  by  McClintock  [5],  and  the 
solution  of  eq .  (24)  for  the  case  b  =  03  .  Except  at  small  volume  fractions  all 
the  curves  predict  ductilities  much  larger  than  those  experimentally  observed. 
Certainly  there  have  been  many  simplifications  made  which  might  cause  the  dis¬ 
crepancy.  Voids  of  all  shapes  are  oriented  in  all  directions  in  the  center  of 
a  tensile  specimen  neck.  The  increase  in  volume  fraction  of  voids  with  applied 
strain  due  to  both  void  growth  and  continued  void  nucleation  must  be 
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DUCTILITY 


FIG.  8  DUCTILITY  VS.  INITIAL  VOLUME  FRACTION  OF 
VOIDS. 
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affected  by  interactions  among  voids.  McClintock  modeled  the  neck  as  composed  of 
cylindrical  voids  oriented  along  the  three  principal  directions  and  entirely 
neglected  interaction.  Here  interaction  is  included  in  an  approximate  sense  but 
only  cylindrical  voids  parallel  to  the  specimen  axis  are  allowed.  A  feeling  for 
the  effect  of  void  orientation  and  interaction  on  ductile  fracture  can  be  gained 
by  comparing  the  various  results.  The  finite  and  infinite  body  analyses,  curves 
I  and  II,  show  the  significant  acceleration  of  void  growth,  and  thus  decrease  in 
ductility,  due  to  interactions  among  voids.  The  two  infinite  body  analyses, 
curves  I  and  III,  imply  that  voids  transverse  to  the  specimen  axis  will  experience 
larger  growth  rates  than  voids  parallel  to  the  axis.  A  plastic  solution  to  the 
former  case  of  unequal  imposed  stresses  transverse  to  the  void  would  be  a  welcome 
tool  in  making  further  studies  (McClintock  [5]  used  an  extrapolation  from  vis¬ 
cous  behavior  to  obtain  curve  III).  The  method  of  approximating  void  interaction 
by  considering  each  void  growing  in  a  finite  body  of  radius  equal  to  one-half  the 
void  spacing  cannot  be  adequately  judged  until  it  is  built  into  a  complete  model 
of  void  growth  and  coalescence .  This  must  allow  different  orientations  of  voids 
and  a  description  of  unstable  flow  between  voids.  Until  then  the  lower  duc¬ 
tility  predicted  by  curve  II  compared  to  curve  I  encourages  further  thought  along 
these  lines. 

From  the  previous  two  sections  we  see  that  the  difference  in  growth  rates 
of  voids  in  perfectly  plastic  and  hardening  materials  depends  directly  upon  the 
triaxiality  of  stress .  For  the  small  triaxialities  realized  in  the  neck  of  a 
tensile  specimen  it  should  then  be  a  valid  approximation  to  consider  the  voids 
as  in  a  perfectly  plastic  material.  Solving  eq.  (7)  for  N  =  .0  gives  a  very 
simple  relation: 
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The  imposed  stress  state  is  given  by  eq.  (22);  the  radius  ratio  rQ/b  is  given 
by  eq.  (25).  Solving  eq.  (26)  for  initial  void  volume  fractions  from  .02  to 
.2  gives  predicted  ductilities  only  a  few  percent  lower  than  those  of  eq.  (24). 
Here  is  a  very  important  simplification  of  at  least  the  initial  growth  process  of 
voids  in  the  neck  of  the  tensile  specimen.  The  initial  growth  can  be  considered 
in  a  perfectly  plastic  material.  Future  research  will  show  if  the  same  simplifi¬ 
cation  is  possible  in  the  entire  growth  and  coalescence  process. 

Discussion 

The  growth  of  single  cylindrical  voids  in  finite  and  infinite  bodies  of 
rigid  power-law  hardening  materials  was  related  to  the  state  of  transverse  stress , 
imposed  strain i  and  hardening  coefficient.  Bounds  to  the  actual  flow  stress 
field  are  presented  which  release  the  void  growth  problem  of  its  dependence 
on  deformation  history.  By  solving  two  extreme  cases,  the  stress  state  necessary 
to  maintain  an  imposed  constant  growth  rate  and  the  growth  rate  resulting  from 
an  imposed  constant  stress  state,  the  behavior  of  voids  in  strain-hardening 
materials  is  seen  to  deviate  from  perfect  plastic  behavior  directly  with  value 
of  transverse  stress  and  hardening  coefficient.  As  imposed  strain  increases  the 
strain-hardening  behavior  asymptotically  approaches  perfect  plasticity  behavior. 
Small  transverse  stress  states  like  those  in  the  tensile  specimen  neck  cause 
little  deceleration  in  growth  rate  from  perfect  plastic  growth  rates. 
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Although  void  growth  is  decelerated  directly  with  hardening  coefficient 
with  stress  state  constant.  Rice  and  Rosengren  [10]  have  shown  that  there  is  a  rise 
in  stress  triaxiality  ahead  of  a  crack  which  is  also  directly  related  to  hardening 
coefficient.  Since  an  increase  in  triaxiality  accelerates  void  growth,  hardening 
could  have  a  net  adverse  effect  on  ductility  for  some  configurations . 

The  theory  is  applied  to  the  fracture  of  a  tensile  specimen  by  considering 
the  interior  of  the  neck  as  composed  of  an  array  of  cylindrical  voids  oriented 
parallel  to  the  axis  with  diameter  to  spacing  ratios  equal  to  the  cube  root  of 
the  initial  volume  fraction  of  voids.  Considering  ductility  to  be  that  strain  at 
which  the  voids'  radii  grow  to  one-half  the  spacing  between  voids,  calculated  as 
if  each  void  were  independently  growing  in  a  cell  of  transverse  dimension  equal 
to  the  void  spacing,  the  resulting  ductility  vs  volume  fraction  curve  was  very 
much  above  the  experimental  curve  of  Edelson  and  Baldwin  [9].  Comparing  results 
with  McClintock  [5]  it  was  concluded  that  voids  normal  to  the  direction  of 
straining  grow  faster  than  voids  parallel  to  the  axis.  Interaction  among  voids 
is  seen  to  be  a  significant  factor  in  void  growth  from  nucleation  to  coalescence. 
The  need  for  a  more  complete  model  of  void  growth  and  coalescence  is  noted. 

Such  a  model  would  allow  random  orientation  and  both  stable  and  unstable  inter¬ 
actions  among  voids.  To  keep  the  analysis  manageable  it  seems  that  only  cylin¬ 
drical  and  spherical  configurations  would  be  included.  The  problem  of  the  growth 
of  a  cylindrical  void  under  unequal  transverse  stress  components  must  be  solved 
if  one  chooses  to  use  the  cylindrical  configuration.  The  results  for  the  growth 
of  a  spherical  void  in  a  perfectly  plastic  material  [7]  might  be  utilized.  The 
finite  body  approximation  to  interaction  could  be  tested  in  such  a  model.  Of 
course  the  instability  of  flow  between  coalescing  voids  is  another  problem  to 
be  solved. 
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The  tensile  specimen  fracture  problem  was  solved  with  the  approximation 
of  perfect  plastic  material  with  near  identical  results  to  the  strain-hardening 
solution.  This  is  a  result  of  the  small  decelerating  effect  of  hardening  at  the 
low  triaxial  stress  states  present  in  the  neck.  This  encourages  future  work 
on  the  tensile  test  to  use  the  perfect  plastic  approximation. 
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